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We study preheating in models where the inflaton has a non-canonical kinetic term, containing 
powers of the usual kinetic energy. The inflaton field oscillating about its potential minimum acts 
as a driving force for particle production through parametric resonance. Non-canonical kinetic 
terms can impose a speed limit on the motion of the inflaton, modifying the oscillating inflaton 
profile. This has two important effects: it turns a smooth sinusoidal profile into a sharp saw-tooth, 
enhancing resonance, and it lengthens the period of oscillations, suppressing resonance. We show 
that the second effect dominates over the first, so that preheating with a non-canonical inflaton field 
is less efficient than with canonical kinetic terms, and that the expansion of the Universe suppresses 
resonance even further. 



I. INTRODUCTION 



After a sufficiently long period of inflation the Uni- 
verse would be cold and devoid of observable matter. 
The energy responsible for driving inflation is trapped 
in the (nearly) homogeneous inflaton field <j>. In order for 
observable matter to emerge from the post-inflationary 
Universe, the inflaton field must couple to additional de- 
grees of freedom in a way that the inflationary energy is 
dumped into observable matter through a process known 
as reheating PHH . For the purpose of reheating, the in- 
flaton can couple to other scalar fields, fermions, or gauge 
fields. 

If the inflaton couples to bosonic fields, such as other 
scalar fields, novel condensation effects can take place. 
In particular, because there is no exclusion principle, the 
inflaton field can transfer a large amount of energy to 
the reheating field x (the "reheaton") in a process that 
is far from equilibrium. Such enhancements result from 
non-linear resonance effects due to the interaction be- 
tween the inflaton and reheaton. For example, for the 
Lagrangian 



C = \{d4) 2 -v(4) 
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the equation of motion for large-scale modes of x ( ne ~ 
glecting the expansion of the Universe) becomes that of 
a harmonic oscillator with a time-dependent frequency, 



X (t) + (k 2 +g 2 cb(t) 2 ) X (t) = 



(2) 



where k is the comoving wavenumber of the X field 
and 4> — <fi(t) is the time-dependent background solu- 
tion for the inflaton. At the end of single field infla- 
tion the inflaton oscillates about the minimum of its po- 
tential, so that the time-dependent frequency u>(t) 2 = 
k 2 + g 2 <fi(t) 2 oscillates with time. It is well-known that 
a harmonic oscillator with an oscillating time-dependent 
frequency can exhibit resonance effects, where the am- 
plitude x(t) °f the oscillator grows exponentially with 
time. This post-inflationary exponential growth of cer- 
tain long-wavelength modes of X is dubbed "preheating" 



and has been studied in a variety of different con- 
texts. See [9H1H for some reviews of this extensive lit- 
erature. Most importantly for this paper, the majority 
of these studies (with the exception of [12|]) have focused 
only on quadratic kinetic terms for the effective theory 
of reheating. 

This approach towards inflationary and reheating 
model building neglects, however, the fact that these sim- 
ple Lagrangians arc really effective field theory (EFT) 
descriptions, only valid at sufficiently low energies. In 
particular, these Lagrangians should be understood as 
having been obtained by integrating out physics above 
some scale A at which new physics (such as new fields 
or new interactions) become important. The effects of 
physics above this energy can be parameterized in the 
EFT through non-renormalizable operators suppressed 
by powers of the scale of new physics, 



Ceff = C 



ECn 
C ™ An-4 



n>4 



(3) 



As an example, consider the two-field Lagrangian, 
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with (j> the inflaton field and p some heavy field with mass 
A. For energies below A, we can integrate out p at the 
classical (tree) level to obtain the effective Lagrangian 
00, 



-e// 
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{dcj>Y 
A 4 



V(<f>). 



(5) 



The effective theory now contains a new contribution to 
the kinetic part of the action for the inflaton. More gener- 
ally, one can consider as a low energy EFT a Lagrangian 
for the inflaton of the form 1 



C e ff - £eff(X,(j>) , 



(6) 



1 This is, of course, not the most general effective theory of the 
form J3j. In particular, we have omitted terms involving higher 
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where X = — i(<9</>) 2 . Inflation with Lagrangians of this 
type can have novel features, such as a speed limit on the 
motion of the homogeneous inflaton and a sound speed 
of perturbations less than one c 2 < 1, that have impor- 
tant implications for inflationary perturbations and mod- 
els [US]. 

It is the post-inflationary dynamics of © that is most 
interesting to us here. In particular, we will couple the 
non-canonical inflaton in ^ to a canonical reheaton field 
through a quartic interaction: 



C pre = C eff (X, <f) + -{dxf - -<?W - ^m 2 x X 2 ■ (7) 



The equation of motion for the reheaton is still of the 
form @. However, since the resonance arising from the 
time-dependent harmonic oscillator ^ is a non-linear ef- 
fect, it is sensitive to the precise profile of the inflaton cj>(t) 
as it oscillates about its potential minimum. The profile 
of the oscillation is in turn sensitive to the non-canonical 
kinetic terms in ([7]). Thus, the modified dynamics from 
non-canonical Lagrangians can play an important role in 
the physics of preheating. 

In this paper, we will discuss the implications of an 
inflaton sector with non-canonical kinetic terms for pre- 
heating. We first discuss in Section [TT| the dynamics of 
a non-canonical inflaton oscillating about its minimum. 
A fairly generic feature is the existence of a speed limit 
on the motion of the inflaton, restricting how fast it can 
move in field space. When the oscillating inflaton satu- 
rates this speed limit, its profile is no longer sinusoidal 
but instead takes a saw-tooth form. In Section IIIII we 
show how the equation of motion for reheaton pertur- 
bations can be recast into a form of Hill's equation, for 
which well-known techniques exist for finding resonance 
bands. Two competing effects lead to modifications of 
the standard theory of preheating: non-canonical kinetic 
terms lead to a "sharper" inflaton profile, enhancing res- 
onance, while the period of oscillation is lengthened due 
to the speed limit, suppressing resonance. The net result 
is that the latter effect dominates, so that the resonance 
for non-canonical kinetic terms is less efficient than its 
canonical counterpart. After illustrating the effects of the 
expansion of the Universe, we summarize our results in 
Section ITVl The Appendix contains further details about 
methods for computing the properties of parametric res- 
onance. 



derivatives. This can be done self-consistently as long as the 
higher derivatives are small for the physics we are interested in, 
which we will argue is the case. See also [T3l for further analysis 
of the validity of the truncation 10 in the context of inflation. 
It is also worthwhile to note that certain Lagrangians of the 
form (O are protected against corrections by powerful non-linear 
symmetries Hal . 



II. NON-CANONICAL KINETIC TERMS 

Let us first focus on the implications that non- 
canonical kinetic terms of the form © have on the mo- 
tion of the inflaton field oscillating about its minimum. 
Inflationary Lagrangians of this form with no potential 
energy have been proposed as an alternative to potential- 
dominated inflation, and are dubbed "k-inflation" 
More generally, however, we expect the inflationary sec- 
tor to have both kinetic and potential energy; in this case, 
the effective Lagrangian C e f / (A, <f>) can be seen as an 
extension of standard slow-roll inflationary models with 
non-canonical kinetic terms. These non-canonical La- 
grangians have a number of interesting properties that 
make them attractive to inflationary model building and 
phenomenology. On the model-building side, the non- 
canonical kinetic terms significantly modify the dynam- 
ics of the inflaton so that the speed of the inflaton can 
remain small even when rolling down a steep potential 
[H|,[l8j]. In many cases, this results in a "speed limit" 
for the inflaton, namely a maximum for the speed of the 
homogeneous mode of the inflaton. This is a great ad- 
vantage to model building, since then it is less necessary 
to fine-tune the inflationary potential to have very flat 
regions supporting slow roll inflation. Inflation with non- 
canonical kinetic terms can also lead to interesting signa- 
tures in the CMB, such as observable non-gaussianitics 

We will choose to work with inflationary Lagrangians 
of the "separable" form (discussed in more detail in [l3], 
HI): 

vn+l 

c eff (x, 4>) = J2 c n ^r - v{<j>) = p(x) - v(cf>) , (8) 

n>0 

where we have written p(X) as a power series expansion 
in Xj A 4 , with A some UV energy scale. Since A is typ- 
ically the mass scale of some heavy sector we have inte- 
grated out, we require it to be larger than the mass of the 
inflaton A ^> in order for such an effective field the- 
ory perspective to make sense. As a further restriction, 
we will take cq = 1, so that for small A/A 4 <C 1, the La- 
grangian reduces to the usual canonical Lagrangian. In 
order for the power series expansion to make sense, the 
series must have some non-zero radius of convergence R 
so that the series converges in the domain of convergence 
A/A 4 e [0, R), with R < 1 (it is not necessary that the 
series itself converge at the boundary). Finally, we will 
also require that the first and second derivatives of the 
power series are positive , J^j- > so as to guaran- 
tee that we satisfy the null energy condition and have 
subluminal propagation of perturbations [IH I22I - I24I ] . 
In a FLRW background 

ds 2 = -dt 2 + a(t) 2 dx 2 (9) 

the scale factor a(t) is driven by the energy density of the 
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homogeneous inflaton <f>(t), 



3M p 2 



2X^-p(X) + V(4>) 



(10) 



where now X = \(j) 2 ■ The equation of motion for (f>(t) in 
this background becomes, 



<3H4>c 2 s 



dV 



d4> dp/dX 



= 



where the sound speed c 2 , defined by 



1 + 2X 



d 2 p/dX 2 
dp/dX 



(11) 



(12) 



is so called because it also appears as the effective speed 
of perturbations of (j> about this background. We will 
ignore the expansion of the Universe for now, dropping 
the Hubble friction term in (JTTJ) , and will return to the 
effects of expansion of the Universe in Section MI CI 

For simplicity, we will take the potential V(<f>) of the 
inflaton to only consist of a mass term V(cf>) = ^rn^cj) 2 . 
This could be the entire inflaton potential, as in chaotic 
inflation, or just the form of the potential near its mini- 
mum. The difference is not particularly important, as 
we are primarily concerned with the phase when the 
homogeneous inflaton is oscillating about its minimum, 
where the quadratic form of the potential will be suffi- 
cient. For a scalar field with a canonical kinetic term in 
this potential, the equation of motion is that of a simple 
harmonic oscillator: the inflaton oscillates sinusoidally 
4>(t) = $sin(m0i). 

The behavior of a non-canonical kinetic term, however, 
can be qualitatively different. As with the canonical case, 
the potential provides a force that accelerates the infla- 



ton. However, now the effective force 



av 



is also 



90 dp/dX 

a function of the speed of the inflaton. Recall that the 
kinetic term p(X) is a series which is only defined within 

a finite radius of convergence, < -R. As the inflaton 
speed approaches the radius of convergence \<f>\ — > V2RA 2 
the series p(X) or its derivatives dp/dX,d 2 p/dX 2 may 
converge or diverge, depending on the precise form of the 
series chosen. Notice that if the second derivative of the 
series d 2 p/dX 2 diverges faster than the first derivative 
dp/dX at the boundary of the domain of convergence, 



so that 



dp/ax 



— > 0, then the effective force 



IX d 2 p/dX 2 

vanishes. Said another way, as the force from the poten- 
tial increases the speed of the inflaton, the non-canonical 
kinetic dynamics modify the effective force felt by the 
inflaton so that the effective force, and thus the acceler- 
ation of the inflaton, vanish as the inflaton approaches 
the radius of convergence. Thus, Lagrangians for which 
d 2 p/dX 2 diverges as <j) approaches the radius of conver- 
gence have a speed limit, such that \<f>\ < 4> m ax = V2RA 2 . 
Certainly, as the inflaton approaches the speed limit, we 
are approaching the boundary of validity for the EFT 




FIG. 1: The motion of the homogeneous inflaton <f>(t) about 
the minimum of the (quadratic) potential differs depending 
on whether the kinetic terms are canonical or non-canonical. 
For canonical kinetic terms the motion is sinusoidal. For non- 
canonical kinetic terms the motion approaches a saw-tooth 
with a much longer period as the inflaton saturates its speed 
limit. This form of the inflaton profile is universal for La- 
grangians that have a speed limit, except for a very small 
region near the turning point which is unimportant for pre- 
heating. 



([HI). Thus, it is advantageous to have a symmetry that 
protects the form of the Lagrangian against further cor- 
rections as this threshold is approached [l5j]. Further, re- 
quiring perturbivity of the inflaton perturbations places 
a bound on the minimal sound speed c 2 , thus restricting 
how close to the boundary of the EFT we can go [25|-[27| . 
Fortunately, there is a window where both the system is 
perturbative and the inflaton is close to its speed limit, 
so that the effects we are interested in below are still 
present. 

For potentials that are sufficiently steep, the inflaton 
quickly attains its speed limit and stays there until it 
reaches the other side of the potential, where it deceler- 
ates and changes direction, as shown in Figure [TJ This 
has two important consequences that will play a very im- 
portant role in preheating: 

(a) The profile becomes sharper, with a smooth sinu- 
soidal profile turning into a saw-tooth profile. 

(b) The speed limit slows the inflaton down, lengthen- 
ing the period. 

These two effects can be clearly seen in Figure [1] As 
mentioned earlier, the form of the potential about the 
minimum is in principle unrelated to the form of the po- 
tential during inflation. This implies that the behavior of 
the inflaton oscillating about its potential minimum can 
be dominated by the non-canonical kinetic terms, inde- 
pendent of whether inflation itself is dominated by these 
terms. 

When the inflaton saturates the speed limit its profile 
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along one of the legs is approximately linear in time 2 : 



(j)(t) « V2i?A 2 (i- t,), 



(13) 



where tj is the time when the inflaton crosses zero. For 
small A this linear approximation of the inflaton is valid 
up to the turning point \<p\ « so the period of oscilla- 
tion is 



4$ 



JVCJi 



'2RK 2 



(14) 



In order for the non-canonical 



with t G h 2 ^, Ii f M ] 
kinetic terms to be important, this must be much larger 
than the canonical period of oscillation Tqr = 2ir/m<fi 
(otherwise the inflaton would not reach the speed limit 
during oscillation), so 



Tncr ^> Tcr 



A « { R 



1/4 



(15) 



This provides a precise condition for the small A limit 
where the system is very non-canonical. Together with 
the requirement that the effective field theory description 
make sense A 3> m^, we have the following regime for A 
where both the EFT description makes sense and the 
speed limit is saturated: 



A 

1 < — < 

rritj, \ R 



1/4 



(16) 



Clearly, for too small of an amplitude $ of the initial 
oscillation this cannot be satisfied. Since the expansion 
of the Universe and the transfer of energy from the in- 
flationary to the reheating sector cause the amplitude of 
oscillation to decrease over time, eventually the inflaton 
no longer saturates the speed limit and the inflaton pro- 
file behavior returns to the canonical limit. 

Perhaps the most well-known non-canonical La- 
grangian that leads to a speed-limit for the inflaton is 
that of the DBI Lagrangian [l8[ (with constant warp fac- 
tor, see [H): 



P(X) 



-A 4 



2X 
A 4 



(17) 



This Lagrangian arises by considering the motion of a 
space-filling D3-brane in a compact space, with <f> taking 
the role of a transverse coordinate of the D3-brane. In 
this case the speed limit has a nice geometrical interpre- 
tation - it is just the effective speed of light for motion 
the extra dimensions (the effective speed of light is not 
necessarily unity if the metric on the internal space has 



non-trivial warping). We can also represent (fT7|) as a 
power series representation in powers of X/A 4 as in © 
with a radius of convergence R — 1/2. But the condition 
for obtaining the speed limit behavior can be satisfied 
by a much larger set of Lagrangians, not just the DBI 
Lagrangian. For example, Lagrangians of the form (ljjj . 



P(X) 



-A 4 



1 XV 



with R < 1 (which includes {jTTJ) for R = 1/2), or 



P(X) 



-A 4 



(18) 



(19) 



are all of the form ((SJ) (with radii of convergence R, 1 
respectively) and all lead to speed-limiting behavior for 
<p. Importantly, though, the details of preheating driven 
by such fields will be insensitive to these different choices 
of Lagrangians - as long as there is a speed limit, the 
profile tfi(t) will be that of the solid line in Figure [TJ 



III. PREHEATING WITH NON-CANONICAL 
INFLATION 

A. Floquet Theory of Resonance 

In the previous section we described the effects of non- 
canonical kinetic terms on the profile 4>{t) of an inflaton 
oscillating about the minimum of its potential, ignoring 
the coupling of the inflaton to the reheaton field. Now, 
we will consider the impact on preheating of the coupling 
between the inflaton and reheaton sectors as in , with 
the inflationary Lagrangian given by ([8]) , 



C pre =p(X) - -m^c 



\(dx) 2 - \m\x 2 ~ \g W(20) 



where we will assume that p(X) gives rise to a speed 
limit. 

While the inflaton is dominated by its spatially homo- 
geneous mode 4> = 4>(t), which is oscillating with period 
T, the reheaton is assumed to have a vanishing back- 
ground vacuum expectation value (VEV) (x) = 0. Thus, 
we will consider fluctuations of the reheaton about the 
vacuum, which may be decomposed into a set of Fourier 
modes 



X = Sx(x, t) = 



d 3 k 
(2tt)3 



+ e-*-5 X fc(t) .(21) 



Ignoring for now the expansion of the Universe, the equa- 
tion of motion for the reheaton fluctuations is 



2 This implies that the acceleration is approximately zero, so 
higher derivative corrections to the effective Lagrangian {(8} are 
very small, as discussed in footnote 1. 



6 X k(t)+(K + g 2 <p(t) 2 )6xk(t)=Q 



(22) 



where K = k 



X' 



This is easily recast into a driven 



harmonic oscillator known as Hill's equation [29J, after 
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a redefinition to a dimensionless time coordinate r = 
Ait (t — t )/T (a prime denotes a derivative with respect 
to r): 



S X k(T)" + [A k + qF(T)]6xk(T)=0, 



(23) 



where F(t) is a 27r-periodic function symmetric about 
r = 0, satisfying f_ n F(t)cLt — 0. Floquet's theorem 
(see [29j]) states that solutions to (f2"B"| have the form 



*X*(r) = e feT 5 (r) + e-^ T 52 (r) 



(24) 



where g(r), <?2(t) are oscillating solutions, and the Flo- 
quet growth exponent fik depends on and q and is 
complex in general. For certain ranges of Ak and q, 
known as resonance bands, the real part of the Floquet 
growth exponent is non-zero leading to exponentially 
growing solutions: 



(25) 



In order to compare models of preheating with differ- 
ent oscillating profiles <p(t) it is more convenient to pa- 
rameterize the growth exponent in terms of physical time 
t: 



(26) 



where (ih = fik-f-i with T the period of oscillation of 
the inflaton. The relation between the two growth ex- 
ponents is that p,k represents the growth per oscilla- 
tion of the inflaton, while fik represents that growth per 
unit of physical time, which takes into account effects on 
the overall growth due to changes in the period. The 
physical growth exponent fik is the appropriate quantity 
to use to evaluate the rate of growth compared to the 
rate of expansion of the Universe. Particle production 
is tracked by the number density of particles, con- 
structed as the energy of a mode divided by the effective 
2 + m 2 + g 2 (j>(t) 2 , and also scales with 



mass m 



eff,k ~ * t >" x 
the Floquet exponent as 



m eff,k ( \8Xk(t)\' 



L eff,k 



(27) 



For any given oscillating profile <p(t) for the inflaton, 
with corresponding F(t) in (1231) . it is possible to numer- 
ically determine the resonance bands; the procedure is 
outlined in Appendix [A] We distinguish two opposite 
regimes: 

1) CR: Canonical Reheating, where the reheaton is 
coupled to an inflaton with a canonical kinetic term 
(or equivalently where A 2 >> <f> at all times so 
that the inflaton effectively behaves canonically) , 



3 Without loss of generality, we will take the real part of the Flo- 
quet exponent to be positive. 



2) NCR: Non-Canonical Reheating, where the re- 
heaton is coupled to an inflaton with a non- 
canonical kinetic term p(X), such that the inflaton 



speed approaches the speed limit 
oscillates about its minimum. 



as it 



First, let us write the equation of motion for the re- 
heaton fluctuations Sxk(t) in the form (J23J) for the CR 
case. As discussed in the previous section, an inflaton os- 
cillating in a quadratic potential with a canonical kinetic 
term has a sinusoidal profile <j>(t) — ^ sm{m^t). The re- 
heaton equation of motion can then be recast into the 
form of a Hill equation ([2"B"]) with the identifications: 



2.T.2 



.4 



2^ + 5 ^$ 



8m 2 ' 
r = 2m*i, F(t) 



qcR 



8m 2 



(28) 
(29) 



Hill's equation in this form is more commonly known as 
the Mathieu equation, and the resonance bands, plotted 
in (K, $) space in Figure [2 take the familiar form from 
previous studies of preheating d, @ . 

In the opposite limit (NCR), when the inflaton has a 
non-canonical kinetic term with a speed limit |<^| ma x = 
%/2~R A 2 , the profile </>(i) becomes a saw-tooth, as in Fig- 
ure Q] Along one of the "legs" of this profile, the inflaton 
is linear in time as in (|13[) , so that the reheaton equation 
of motion becomes: 



Sj(k(t) + {K + 2Rg 2 A 4 t 2 ) 5 X k(t) = 



(30) 



This can also be rewritten in the form of a Hill's equation 
(T23l) with the identifications: 



A,= 



K<S> 2 



1 <? 2 $ 4 



g 2 & 



2i?7T 2 A 4 ' 3 2i?7T 2 A 4 ' QNCR 2i?7T 2 A 4 ' 



2RA 2 tt 



(31) 



The resonance bands for this form of Hill's equation are 
shown in Figure [5J The two regimes can easily be con- 
nected numerically. 

There are two main regimes of interest of (|2"3")1 for the 
resonance bands: narrow resonance q 1 and broad res- 
onance q 3> 1. For g 2 $ 2 /r7i 2 i <C 1, the CR scenario is in 
the narrow resonance regime. When non-canonical ki- 
netic terms are important, however, the parameter qNCR 
in the corresponding Hill's equation (|31[) is enhanced rel- 
ative to the canonical case 



qNCR 



9 



2Rtt 2 



TO0 



2 m 2 $ 2 



4 77J 2 $ 2 

Rn 2 A 4 



qcR 



(32) 



The enhancement factor m 2 i <i> 2 /A 4 3> 1 is large in or- 
der for the inflaton to saturate the speed limit, so unless 
qcR is correspondingly small, when the CR scenario is 
in the narrow resonance regime the NCR scenario is in 
broad resonance. This has important physical implica- 
tions, since then not only is the growth per period fik 
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FIG. 2: Resonance bands for reheaton perturbations as a function of the scale K and initial amplitude of the oscillating inflaton 
$ with g — 0.005 in the case when the reheaton is coupled to an inflaton with canonical (CR scenario, left) and non-canonical 
(NCR scenario, right) kinetic terms. 



larger for non-canonical kinetic terms, but there is also 
growth over a much larger range of scales. We will see 
in the next section, however, that this enhancement is 
overwhelmed by suppression of particle production due 
to the lengthening of the period. 

The resonance due to (1221) depends on the physical 
wavenumber of the fluctuation k; most particle produc- 
tion occurs when the effective mass K + g 2 <fi 2 (t) van- 
ishes. This implies that the resonance is most efficient 
at large scales K ~ 0, as can be seen in Figure [2] In 
practice, however, we cannot work on scales larger than 
the Hubble radius while neglecting metric perturbations, 
so k > H. In addition, a non-zero bare mass m x for 
the reheaton also keeps the effective mass from vanish- 
ing. But as long as we work on sufficiently large scales 
(and with a sufficiently large initial amplitude) so that 
K -C g 2 & 2 , the maximum of resonance at K — will 
be a good approximation for the maximum resonance at 
large scales. 

Figure |3] displays the growth exponent fik as a func- 
tion of the initial amplitude $ of the inflaton for large 
scales (namely K — 0). Several features are evident: for 
decreasing A, more resonance bands become accessible, 
because the system enters the broad resonance regime; 
the maxiumum size of the growth exponent in the first 
resonance band decreases as A decreases, reflecting the 
lengthening of the period; and the heights of the reso- 
nance bands for small A decrease with increasing $, also 
due to the lengthening of the period for large initial am- 
plitudes (note that the period for canonical kinetic terms 
is independent of the amplitude) . 

This effect is not limited to the DBI case: Figure 0] 
compares the Floquet exponent for several Lagrangians 
of the forms (Tl8|) and (fT9|) . While the behaviors differ 
slightly as one moves away from the canonical case, they 
converge again in the regime where inflaton oscillations 
saturate the speed limit. The black dotted line in Figure 
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FIG. 3: Floquet growth exponent fj,o per unit of physical time 
for K = for three different values of A in DBI inflaton- 
driven preheating. Lowering the speed limit A 2 of the inflaton 
greatly reduces the strength of parametric resonance because 
the period of oscillation (I14|) increases as 1/A 2 . The addi- 
tional suppression of the resonance due to the dependence of 
the period on the amplitude of inflaton oscillations $ is also 
evident. 



[4] shows the Floquet exponent when the period lengthen- 
ing is not taken into account: as the inflaton trajectory 
approaches the saw-tooth shape, energy is injected into 
the reheaton field over a longer amount of time, allow- 
ing for more particle production. This saturates when 
the slope <f> approaches the constant speed limit for most 
of the period. This enhancement is clearly subdominant, 
however, when compared with the suppression of particle 
production from the elongation of the period itself. 

We close this section by comparing our results with 
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previous explorations in the literature. Ref. [12| found 
an expression for small </>/A 2 for the growth exponent fik 
in the case of a DBI inflaton by perturbing the canonical 
equation of motion. Up to order <fi/A 2 and for K = 0, 
this corresponds to: 




1/2 
; 



(33) 



with 



2m 2 



1 



32A 4 



^2 

4m 2 



1 



3$ 2 m 2 



8A 4 / 

(34) 

This result is illustrated in the right panel of Fig. SJ For 
/A < 1 it provides a good approximation of the gain in 
particle production one would expect from DBI inflation 
were it not for the lengthening of the period, which the 
authors of [l2| correctly identified as "DBI friction" but 
did not quantify. We are interested in the opposite limit, 
however, where the non-canonical kinetic terms are more 
than just a small perturbation. 



B. Non-Canonical Limit 

In the non-canonical limit, the equation of motion for 
the reheaton perturbations Q30p can be written in a sim- 
ple and suggestive form by making a redefinition to a 
different dimensionless time coordinate than we consid- 
ered earlier: 



r = (g V2^A 2 ) V2 t = (2Rg 2 ) 1/4 A t . 
The reheaton equation of motion then becomes 
d 2 Xk(r) 



where k 2 



dr 2 

k 2 +m [ 



/2ft S A 2 



+ W + t 2 ] x/t(r)=0, 



(35) 



(36) 



The time range of the new time 



variable is < r < 4f where At = T A (2g 2 R) 1 / i A = 

(8g*/R)V*($/A). 

The benefit of making this redefinition is that the prob- 
lem of broad resonance can be mapped to the problem 
of scattering of a particle with energy k 2 off of an in- 
verted parabolic potential with an effective Schrddinger 
equation [!, H|: 



d 2 Xk(r) 
dr 2 



+ ( K 2 -V eff (T))xk(r)=0. 



(37) 



This scattering problem can be solved using standard 
techniques (see Appendix [Bjl . leading to the result: 



2tt 



1 



2 . - .„ cos 6 tot 



(38) 



in terms of the reflection and transmission amplitudes 
Rk,Dk and the phase of the reheaton mode 8 to t- For an 



inverted parabolic potential, the reflection and transmis- 
sion coefficients are known, so the average growth index 
simplifies to be 



/2 fc = -!- In 

Z7T 



1 + 26-™' -2sin(9 fot e 



(39) 



For large scales k 2 ps this has a maximum/average value 
of ji k re (0.28,0.175) 0. 

The average growth index Jx^ ps 0.175 only gives the 
typical growth over one period of oscillation for broad res- 
onance, and is the same regardless of whether the inflaton 
dynamics are canonical or non-canonical. In particular, it 
does not depend on the type of non-canonical Lagrangian 
used - any non-canonical Lagrangian that leads to a speed 
limit for the oscillating inflaton will have the same aver- 
age growth factor over one period flk- The reason the 
average growth index does not depend on whether the 
system is canonical or non-canonical is because we are 
working in the broad resonance limit, where only the 
slope of the inflaton profile (in dimensionless coordinates) 
at the point where the inflaton crosses zero matters. 

However, as discussed before, the true comparison of 
the particle production between the canonical and non- 
canonical models should be the growth /i^ over some fixed 
physical time period At 



n/c cx e 



(40) 



where T is the period, so that the true growth index is 
the average growth index over one period divided by the 
length of the period, \ik — jsk^/T. This leads to the 
main result of this section: since the period for an os- 
cillating inflaton with non-canonical dynamics is much 
longer than the period for an oscillating inflation with 
canonical dynamics T/vcft 3> Tcr, we have much more 
growth in a canonical system over some fixed physical 
time At, (Mfe)cft (Mfe)jvCft- ^ n particular, the ratio 
between the two should fall off as the ratio of their peri- 
ods, 



(Mfc)jvcft _ 



Cft 




2 <I>r7i0 



«1, 



(41) 



for decreasing A at k — 0. Indeed, this is precisely what 
is found using numerical techniques, as shown in Figure 
® 

Finally, we comment on the possiblity for non- 
canonical kinetic terms to enhance resonance by turning 
a putative CR system in narrow resonance into a NCR 
system with broad resonance. A CR system in narrow 
resonance with geft < 1 has (fik)cR ~ qcR"m<p @- As 
discussed in Section IIII A[ the same system in NCR has 
qNCR ~ qcR.f 2 , with / = $m /(V^RA 2 ) > 1. But as 
just discussed the growth exponent in broad resonance 
for NCR is a fixed number divided by this enhancement 
factor, (fik)NCR. ~ O.lm^//. Thus, in order for the 
NCR system to be in broad resonance qNCR <^ 1, we have 

1/2 

that {Hk)NCR S5 O.lm^gp^j, representing an enhancement 
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FIG. 4: Floquet exponent fj, of the first resonance band for K = for g — 0.005 (left) and g = 0.5 (right). A >> corresponds 
to the canonical case CR. Parameter R in the left-hand panel refers to (JTHJ) , and log refers to ()19|) . In both panels the black 
dashed line is the growth of particle production per period, whereas the blue solid line is the true growth of particle production 
per physical unit of time. Additionally, in the left panel the magenta dotted line is the analytical result (|41[) . showing good 
agreement with the numerical results in the regime where the speed limit is important. In the right panel we have also included 
the approximate result from [121 ] (dotted red line, (|33p ). which did not include the suppression from the speed limit. As 
discussed in Section UT] the effective theory breaks down for > A and results should therefore not be trusted beyond this. 



from CR narrow resonance by a factor of 0.1/ qJ R - While 
this enhancement can be large if the original CR system 
is very far in the narrow resonance regime, the growth 
factor is still quite suppressed. 



C. Expansion of the Universe 

Up to now we have ignored the expansion of the Uni- 
verse. Including expansion has several effects onpreheat- 
ing which have been well-studied for CR [!, H, Hjj : 

i) The amplitude of inflaton oscillations decreases due 
to Hubble damping inversely proportional to time 
$(t)~f-i. 

ii) The parameters of Hill's equation (Ak,q) become 
time dependent. In particular, qcR ~ t due to 
the time dependence of the inflaton amplitude. 

hi) The physical wavenumber k p h ys = k/a(t) redshifts 
with the expansion, so each comoving mode only 
spends a finite amount of time in a resonance band. 

iv) The resonance is stochastic [3j]. 

v) Preheating is ineffective when the rate of parti- 
cle production drops below the expansion rate, 
e.g. when fik < H. For narrow resonance this con- 
dition cannot be satisfied, thus particle production 
can only occur in the broad resonance regime. 

We will now make some simple estimates of how these 
effects will change for NCR. 



First, the Hubble friction term may be removed from 
the inflaton equation of motion (fTTj) by rescaling the 
inflaton field by powers of the scale factor <f>(t) = 
$o</>(*)M*) 3/2c ". so that 4>(t) is a fixed unit-amplitude 
oscillating function, and $(f) = $o/ a (X) 3/ ' 2Cs is the time- 
dependent amplitude. When the period of oscillation of 
the inflaton is much smaller than the timescale for the 
expansion of the Universe, the energy density can be av- 
eraged over several oscillations ~ rn^cj) 2 ~ et(i)~ 3 / c =. 

The scale factor thus grows as a(t) ~ t 2 / 3Cs , so that the 
amplitude falls is inversely proportional to time ~ 
*f>o/t as with a canonical kinetic term. The decaying in- 
flaton profile in an expanding background is shown in 
Figure [5] It is worth noting that in the far non-canonical 
limit the period of the inflaton can be larger than the 
Hubble time due to the period lengthening effect. When 
this occurs, Hubble friction damps a significant amount 
of the amplitude in the first oscillation, quickly bringing 
the system back to the previous case. 

For a fixed cutoff scale A, the condition for the os- 
cillating inflaton to saturate the speed limit (fPo)) is vio- 
lated once the amplitude of oscillation drops below some 
critical amplitude. The number of oscillations of the 
inflaton scales with time N ~ t/T, so that the time- 
dependent amplitude drops below this critical amplitude 
after N ~ m^&o/A 2 oscillations, where $o is the ini- 
tial amplitude. Thus the inflaton no longer saturates the 
speed limit after N ~ m^&o/A 2 oscillations, and behaves 
canonically after this time. 

Another important effect is the dependence of the 
Hill's equation parameters (|31l) on time. Since Qncr de- 
pends on the fourth power of the amplitude, it quickly 
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-0.01 - 
-0.02 - 




FIG. 5: Expansion of the Universe causes the amplitude of the 
inflaton to decay inversely proportional to time, as shown here 
for A = 10 m<$> and — 10~ 5 M p with an initial amplitude 
of $o = 0.05M P . The dashed line indicates the t _1 behavior 
of the decaying amplitude. 

becomes small at a much faster rate Qncr ~ t~ 4 than 
it would in CR, qcR ~ t . In particular, after only 
N ~ g 1 / 2 $>o/A oscillations, the narrow resonance regime 
where Qncr ^ 1 is attained. As noted above, narrow 
resonance with CR is inefficient for an expanding back- 
ground. Since we argued in the previous section that par- 
ticle production for NCR is suppressed compared to CR, 
narrow resonance for NCR is also inefficient in an expand- 
ing background. This means that after N ~ g 1 / 2 G>o/A 
oscillations the resonance shuts off due to expansion of 
the Universe. The other effects of an expanding Universe 
listed above for CR (iii,iv) do not appear to be qualita- 
tively different for NCR, so we will not comment further 
on them. 

In summary, including expansion of the Universe for 
NCR leads to a decaying amplitude of inflaton oscilla- 
tions inversely proportional to time ~ t _1 . After 
many oscillations the amplitude decreases so much that 
either: (a) the inflaton ceases to saturate the speed limit 
during oscillation, effectively becoming canonical; or (b) 
the Hill's equation parameter qNCR(t) becomes so small 
that particle production cannot compete with the expan- 
sion of the Universe, shutting off the resonance. Which 
outcome dominates depends on the relative magnitudes 
of 771,0 /A and g 1 ^ 2 ] when the former dominates, outcome 
(a) occurs first, and vice-versa. 



IV. CONCLUSION 



some energy scale A > m^. Effective theories of this 
type can give rise to a speed limit for the motion of the 
inflaton cj), as in DBI inflation [l8[. In addition to having 
important implications for inflationary model building, 
the speed limit plays an important role in modifying the 
nature of preheating in the post-inflationary Universe. 
In particular, as the non-canonical terms become impor- 
tant, non-canonical preheating departs significantly from 
the canonical case via three main effects: 

i) The sinusoidal inflaton profile becomes a saw-tooth, 
elongating the fraction of the inflaton period in 
which significant particle production may occur, 
and moves the system from narrow to broad res- 
onance. 

ii) Effect i) is offset by an elongation of the inflaton os- 
cillation period by a factor / = \J^-^r 3> 1. This 
suppresses the amount of \ particle production per 
unit time by 1//. 

iii) Effect ii) affects the competition between \ P r0 ~ 
duction and Hubble expansion, making preheating 
even less efficient in an expanding Universe. 

In general, then, preheating when the inflation has non- 
canonical kinetic terms is less efficient than with a purely 
canonical kinetic term. This implies that if preheating is 
to be important at all in the early Universe, then the 
UV scale of new physics that couples kinetically to the 
inflaton must be sufficiently high that the effective non- 
canonical kinetic terms are negligible. 

We have only focused on non-canonical kinetic terms 
for the inflaton sector in this paper. Certainly, how- 
ever, due to the non-linear nature of the parametric reso- 
nance, it would be interesting to study how non-canonical 
kinetic terms for the reheating sector would affect the 
physics of preheating. We leave this to future work. 
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Preheating in the post-inflationary Universe is the ex- 
plosive production of particles far from thermal equilib- 
rium that occurs as the inflaton field oscillates about its 
potential minimum. We have examined preheating for an 
inflaton sector that has non-canonical kinetic terms aris- 
ing as an effective theory C e ff(X,<f>), for X — i(<90) 2 , 
which may arise from the existence of new physics at 



Appendix A: Computing the Floquet exponent 

There exists a straightforward method of computing 
the Floquet exponent fik as a function of the parame- 
ters Ak and q, which has been extensively covered in the 
literature (e.g. (UMI). 
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We start with Hill's equation of the form (|23l) . The 
periodic function F(r) may be decomposed: 



M 



f(t)= d ^ mT > 



n=-M 



with Ai- defined such that: 



[ F(r)dT = 0. 

J — IT 



(Al) 



(A2) 



Floquet's theorem (see [29]) states that solutions to ([23 
are of the form, 



Xk(T) = e^ T g(T)+e-^ T g 2 (r) 



(A3) 



where g(r) and gi (r) are periodic functions with period 
T, and fik, called the Floquet exponent or characteristic 
exponent, is complex. Clearly when fl k has a non-zero 
real part we have exponential growth of Xk — this is the 
parametric resonance effect. Without loss of generality, 
we will take the real part of jl k to be positive, and so 
we will drop 32 ( T ) since its coefficient is exponentially 
decreasing. 

In order to find solutions of the form (|A3[) . we first 
Fourier expand g(r) 



Xk 



(r)= J2 c„e fe+m)T 



(A4) 



n— — 00 



and plug this Fourier series back into the equation of mo- 
tion ([23| to derive a recursion relation for the coefficients 
c„ in terms of fi k ,A k ,q, d n : 

c n + -7TT — . , 2 . - =0 Vne (-00, 00 . A5 
{{li k + m) z + Ak) 

These recursion relations define a matrix problem 

( C - n \ 

c-i 

B(p k ,A k ,q,d m ) c =0 (A6) 

Cl 
V Cn I 

with n — ¥ 00, where the elements B rs of the (infinite) 
matrix B are given by: 



J 1 if r = s 

B r s — \ qd r - s 1 

I (fik+iry 2 +A k r T 1 S 



(A7) 



The matrix problem (|A6j) implies that the (infinite) 
matrix B is singular, so it must have vanishing determi- 
nant 

A(fL k ,A k ,q,d m ) = \B rs \ = 0. (A8) 
The vanishing of the determinant (|A8|) defines an implicit 
function fb k — p, k (A k ,q,d m ). In practice, we perform a 



Fourier transform of the inflaton profile written in terms 
of F(t) numerically, and evaluate the MxM determinant 
(|A8j) : a matrix size M ~ 100 is sufficient for convergence 
of the first tens of resonance bands. 

Let us briefly comment on the differences in solving for 
the growth exponent p, k using this method for canonical 
and non-canonical kinetic terms. Clearly, F(t) for non- 
canonical kinetic terms pip expanded in Fourier modes 
involves many terms, as opposed to the single Fourier 
mode for canonical kinetic terms. Thus, one difference 
between resonance with canonical and non-canonical ki- 
netic terms is the inclusion of more terms in the Fourier 
expansion of the profile. Naively, then, there are more 
modes in the driving force with which the Xk fields may 
resonate. This is the feature arising from the "sharpen- 
ing" of the profile. However, the other important effect 
of the non-canonical kinetic terms is the extreme length- 
ening of the period of oscillation, thus suppressing the 
resonance in physical time. 



Appendix B: Scattering from a Parabolic Potential 

In many cases of interest, the equation for perturba- 
tions of the reheaton Xk- 



Xk{t)+[k 2 + m 2 x +g 2 ^] X k(t)=0, 



(Bl) 



can be converted into a Schrodinger scattering problem. 
In this appendix, we will outline the equivalence between 
these two perspectives; see @ for more on the correspon- 
dence. We will imagine that <p(t) has a large slope in 
some region for some finite amount of time, as shown in 
Figure [6] Linearizing about this point, we have 



<j> w 0o + A(i - U) 



(B2) 



This linear approximation is valid roughly for some time 
interval 2At, i.e. for tj — At < t < ti + At. Outside this 
interval, the field x must behave adiabatically: its poten- 
tial is slowly- varying and can therefore be approximated 
by a WKB approach. The equation of motion around the 
point of adiabaticity violation becomes: 



Xk(t) + [k 2 + m 2 x + g 2 (fo + X(t - U)) 2 } X k{t) = .(B3) 
We can map this problem into a Schrodinger-like equa- 

1 /2 

tion by redefining the time variable r' = (j) (<f>o + 
\(t — t{)), so that the equation of motion is written: 



d 2 Xk(r') 
dT n 



+ [n 2 +r' 2 ]xk(r / )=0. 



(B4) 



We have defined 



k 2 +ml 



The time range of 
the new variable is t[ — At' < t' < t[ + At', where 
t[ = <j) {g/\) 1/2 and At' = (Xg)^ 2 At. The transformed 
equation (|B4p resembles a Schrodinger equation with the 
effective potential: 



11 



\t[ - At 1 ) 2 = -V = const r' < r/ - At' (Region I) 
V e// = { -t' 2 t[ - At' <t' <t![ + At' (Region II) (B5) 

-(t- + At') 2 = -Vi = const t' > t[ + At' (Region III) 



The setup of this effective scattering problem is shown in 
Figure |6] 

The solution in each region is therefore known: 



for At = T over one period. Thus, the growth exponent 
is, 



A ie -ikir' + B 1 e +lklT ' 
V/, = \ A 2 D V ((1 + i)r') + B 2 Dij ((- 
A 3 e- lk2T ' + B 3 e+ lk2T ' 



(Region I) 
i)r') (Region II)(B6 1 
(Region III) 



1 

2^ 



log 



l + |J?fc 
\D k \ 2 



\R k 



D k 



■ cos Q t 



(B12) 



where fci = y«? 



V 2 



a/k 2 + V£, v = - k 2 ), 



and the functions D v (x) are parabolic cylinder functions. 

In the language of scattering matrices, the ingo- 
ing/outgoing waves are parameterized as 



j+i - 

°i e 



Pi 




(B7) 



where D kl R k are the transmission and reflection coeffi- 
cients, ihea k ,/3 k are normalized such that |a/c| 2 — \/3 k \ 2 = 
1 with the "particle number" defined as n k — \/3 k \ 2 , and 

6{ = J Q T ^k 2 + m 2 , + g 2 (/) 2 dt is the accumulated WKB 

phase up to the time of scattering. Apart from an overall 
normalization, it should be clear that the a? 's correspond 



to the A,-'s and the 



P's to the B/s above. 



From (IB7I) we can write 



-2i8] 



Rk 



(B8) 



which we can use to find the final number of particles 
n^, +1 in terms of the original number of particles n\ (to- 
gether with the normalization condition on the a): 



i+i 

T 



Rk 



Dl 



+2- 



\Rk 



Di 



1 

~\D k \ 2 
■ cos 9tot 



+ n i) n i 



(B9) 



where 9 to t = arga^, — arg/3^ — aigR k — 29 J k is the total 
accumulated phase. In the large nL limit, we can drop 
the first term and simplify the last term, so that we have 



Averaging over many possible initial phases, we can 
consider 9 to t as a random variable; thus, we will take 
cos Qtot ~ on average. 

In order to compute the scattering amplitudes R k , D k , 
one needs to solve the continuity equations across the 
boundaries of Regions I, II and III for the wavefunctions 
(|B6|) . In general, it is not possible to solve these equations 
analytically. However, for special cases the equations can 
in fact be solved. 

The first limit we will take is \t[\ 3> 1; in the orig- 
inal setup, this corresponds to taking 4>o\/g/V^ to be 
large, e.g. the region of linearization takes place very far 
from tj> = 0. In this limit the effective potential becomes 
essentially a step function: 



Veff = 



-V 
-Vi 



(Region I) 
(Region II) 



(B13) 



where \Vq — Vi\ = y/g/X\A<f>\, with Acj) the change in the 
inflaton over the time period At, and \ Vq\ ~ | Vx ] 3> 1- In 
order for the step function approximation to be a valid 
approximation to the inverted quadratic effective poten- 
tial, we need that the change over the step function be 
small, e.g. I^ "^ 1 ! *C 1. The transmission and reflection 



Vi 



coefficients can easily be computed 
A(k 2 + V )(k 2 



Ri 



Vi). 



(2k 2 - 
(V 



Vo + V x ) 2 
-V) 2 



(2 K 2 + V + V 1 ) 2 ' 
so that the growth exponent becomes 

(Vo-V,) 2 



Hk = 



2vr 



1 



V(k 2 + Vo){n 2 + V x ) 



(B14) 



(B15) 



J+ 1 



i + \Rk\ 2 , _ \R k \ „ 

— 1" 2 .„ ., COS 9 tot 



(BIO) 



The average growth exponent over one period is defined 
as 



(Bll) 



With the requirements above that the size of the step be 
small, we have D k 1 and Rk ~ 0, so that ftk ~ m (l) ~ 
0. This is expected - if the inflation has a very large offset, 
the effective mass of the reheaton m 2 ff — k 2 +m 2 x + g 2 (f) 2 ) 
is large, so there will be very little particle production. 

Alternatively, in the limit t[ = 0, the scattering poten- 
tial becomes symmetric about t = 0. For At ^> 1 the 
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K 



v., 



t,+At 



t-At 



FIG. 6: Left: Linearizing <f>(t) w <f>o + A(t — U) about some point (U, 0o). Right: The associated effective scattering problem. 



scattering amplitudes are known 



The implications of scattering in this case is considered 
in more detail in the main text. 



D k = 

Rk = 

where the angle ipk is 



Vl + e"™ 2 
Vl + e™ 2 ' 



,l + i K 2 \ k 2 ( 2 
^ = axgr(^- +y 1 + ln-, 



(B16) 
(B17) 

(B18) 
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